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This is an expository paper presenting the following: (1; the origin, and (2) the ])ro])- 
erties of the Rayleigh distribution; (3) the most efficient estimators of its parameters; (4; 
a test of the hypothesis that a set of observations is from a Rayleigh distribution; (5) the 
distribution of the ratio of two independent Rayleigh variates; and (6) the Rayleigh 
process derived from a normal process. 



1. Introduction 

In son).c problems of physics and ongiiioeriii^'', 
such as the diffusion process, or the distribution of 
the power (or ainplitude) of electromagnetic waves 
received through a scattering niecUuiu, we are con- 
cerned witli the resultant of many (two- or tlu'oe- 
dimensional) raudoiu vectors. Such problems are 
treated in jorobability theory under tlie lieading of 
^^random walk'^ in one, two, three, or, in general, k 
cUmensions. Here, we will confine our attention to 
two-dimensional random walk problems which 
asymptotically give rise to the so-called Rayleigh 
distribution. Unfortunately, in radiowave propaga- 
tion literatin^e the name '^Rayleigh distribution^' is 
used indiscriminately both for the distrit)ution of 
power and of am])litude. To distinguish between 
the two, we will call them. Rayleigh power distribution 
and Rayleigh amplitude distribution, respectively. 
In statistical literature the distinction is clear: the 
Rayleigh power is a multiple of x^ (with 2 degrees 
of freedom), whereas the Rayleigh amplitude is 
a multiple of x- Since the properties of tlie Ray- 
leigh power distribution are nmch neater and better 
suited for statistical estimation and for testing of 
hypotheses than the properties of the Rayleigh 
amplitude distribution, we will study only the 
Rayleigh power distribution. Important charac- 
teristics of the amplitude distribution will be ob- 
tained from the relationship amplitude = (power )^''^ 
The advantages of power over amplitude are roughlv 
tlie same as those of variance over standard devia- 
tion, or of spectral density over its square root. 
For example, the powers are additive and amplitudes 
are not; the smn of two independent x"^ variates is 
again a x" variate witli added degrees of freedom, 
but tlie sum of two x variates is not a x variate. It 
is, tluM'c^fore, suggested that the nam.e ' 'Rayleigh 
distribution, '^ if used without any other qualification, 
be reserved for Rayleigh power distribution, and 
'^Rayleigh variate'^ for a multiple of x^ Avith 2 
degrees of freedom. 



2. Notation 

Orchnarily the same letter will be used to designate 
both a random variable and the argument of its 
probability density or distribution function. Thus 
we may say : a^ is a random varial)le witli i)robability 
density function 2^ (x) and distribution function P(a:). 
If it becomes necessai'y to (hstinguisli b(^tween the 
two to avoid confusion or ambiguity, ('a])ital letters 
will be used for the random variables and h)wer case 
letters for the real variables. Thus we may say : X 
is a random variable and rr{x<X<x-\-dx)=p(x)dx. 

The function P, with one or more arguments, will 
be used generically for any distribution function, and 
p for any density function. Thus P{x) and Piy), in 
general, will not t)e tlie sanu^ functions. The former 
denotes the distribution function of x and the latter 
that of 7/, and the two may be different. 

If x and y are random variables. Ex will stand for 
the expected value of x, and E(x\y) for the condi- 
tional expected value of x given the value of y. Thus 

J 00 /» CO /^ c» 

xp(x)dx= xp(;x,y)dydXj 

— 00 »/— 00»/ — 00 

It is easy to see that, if Ex<C^ , 

/• 00 /^ 00 

= xp{x,y)dy dx=Ex. 

%) — CO^ — CO 

The characteristics of a theoretical distribution, 
i.e., its population parameters, will be designated by 
Greek letters such as /z, a, a. Their estimates may 
be denoted by corresponding Latin letters, with or 
without subscripts or primes. Thus if n denotes the 
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mean (which is also the median) of a normal distri- 
bution, m may denote the sample mean and rrii the 
sample median. 

3. Asjrmptotic Distribution of the Resultant 

Let ri=(ri cos a^, r^- sin at), i=l, 2, . . ., n he in- 
dependent random coplanar vectors where all values 
of ai(i=l, 2, . . ., n) between and 2t are equally 
likely. Let Ti be a series of positive constants, or 
else be positive random variables distributed inde- 
pendently of ai. Then 



Evi cos ai=EriE cos ai=Er 



'2tJo 



cos aidat^O, 



(j\=Er\ cos^ ai=E7JE cos^ cii=-^ Er\. 

Li 



(3.1) 
(3.2) 



Similarly, 



Evi sin ai=0, Er\ sin^ «z=o ^A) Et\ sin a^ cos a^^O, 

(3.3) 

SO that the components of a vector are uncorr elated. 
This does not m.ean, however, that they are statis- 
tically independent. We note that if r^is a constant 
ETi=T,, Erl=rl 

The resultant vector r = {x,y) is given by 



3^=S ^i <^os Gi, 2/=2D ^i sin di. (3.4) 

i=l i=l 

We immediately have 

Ex^Ey^O, a'=^Ex'=Ey'^f]<Tl (3.5) 

i=l 

Applying the Central Limit Theorem for vector 
variates [Cramer, 1951, pp. 215-216, 285-286], we 
obtain the following theorem. 

Theorem. Suppose that p\=Er\<i^ for alii, and 
write 

P — Pi-r • • '-rPn- 



Ij the condition 



lim p/(r=0 



is satisefid, then, noting that x and y are uncorrelated, 
we obtain the asymptotic distribution of (x,y) as 



dP(x,y)- 



2'jra' 



e ^'"^ dx dy, — oo < x, y< oo . 



Thus xjd and yjc are asymptotically independent 
standard normal variates. 

We note that if r^ have identical distributions and 
£rf < 00 ^ or if rj = r, a constant for all i, p/(T=0(n~'^^^)^ 
which tends to zero as n-^ oo . 



If in the equation above we make the transforma- 
tions 

x=r cos a, y=r sin a, (3.6) 

the Jacobian of the transformation is r, and the 
distribution of (r,a) is given by 

dP(r,a)=-,e ^'''dr • ^; 0<r< oo , 0<a<27r. (3.7) 

Hence the resultant phase, a, and the resultant 
amplitude, r, are independently distributed. We 
will write y = 2a^, so that the distribution of r is 



r2 

dP{r)=—e ''dr, 0<r<oo, 
7 



(3.8) 



This distribution is known as Eayleigh amplitude 
distribution. Making further the transformation 



2=zr^ 



(3.9) 



we obtain the distribution of z, called the Eayleigh 
power distribution, as 

dP{2)=^^ e~'^dz, O<0<oo. (3.10) 

7 

We note that 2=x^-^y^is the sum. of two independent 
normal variates, x and y, each with mean zero and 
variance 7/2. Thus 2z/y is a x^ variate with 2 
degrees of freedom [Cramer 1951, p. 233]. 

In case all ri = l in (3.4), Kluyver [1906] showed 
that 

P{r) = rjjMt)rMrt)dt. (3.11) 

If we make the transformation 

t=2n-'^'y 
and observe that 

lim[Jo(27i-i/2y)]«=6-^', 
an asymptotic approximation to (3.11) is obtained as 



P(r)^l-e ^. 

This formula is originally due to Rayleigh [1919]. 
In general, if we start with k dimensional vectors 

-> -^ '^ -> 

Vi^ixn, . . ., Xi^), and r=2]^i = fe, • • -, ^J, then, 

i=l 

under the conditions of the Central Limit 
Theorem, Xi, . . ., a:;^ are asymptotically independent 
normal variates, each with mean zero, and the 
same variance, 7/2. Let z^x\-\- . . . +x\. Then 
2^/7 is a x^ variate with k degrees of freedom 
[Cramer 1951, p. 233]. Thus the asymptotic 
distribution of z is given by 



168 



4. Rayleigh Power Distribution 
4.1. Distribution Characteristics 

Til tliis section we will study some population 
characteristics of the Rayleigh power distribution. 
In the following sections the problems of estimation 
and of testing of hypotheses concerning this distri- 
bution will be considered. 

The probability density function of z is 

2){z)=-e~^^'^, if 2>0; 0, otherwise, (4.1) 
y 

which has a maxhn.um at z = and decreases mono- 
ton ically as z increases. 

The diMribution functioti, Pr(Z<z), is givcMi by 



F(3). 



7 Jo 



ydx=l — e y* 



(4.2) 



If g is a number between and 1, tin* lOOq per- 
centile point, ^q, is defined by the ec) nation 

which has the solution 

f, = -7ln (1-r/). (43) 

In parti(ndar the median 

f 5=.^ In 2=0.69315 7. 

The values of ^Jy for (^ = 0.01, 0.05, 0.10, 0.25, 0.50, 
0.75, 0.90, 0.95, 0.99, are given in table 1. 







Table 1. 


.tq/y for se 


lected 


q 






Q 


0.01 


0.05 


0.10 


0.25 


0.50 


0.75 


0.90 


0.95 


0.99 


Uh 


0.01005 


0. 05129 


0. 10536 


0. 28768 


0. 69315 


1.38629 


2. 30259 


2. 99573 


4. 60517 



If k is any real number greater than — 1 

^if,{z)=Ez''=~ f z'e'ydz^y'^rik+l). (4.4) 
7 Jo 

Thus the mean, ^i(z), and the variance, (t^(z), are 

}xi(z)==y, a\z)=fi2(z)-fii{z)=y\ (4.5) 

Also the mean and variance of r=z'^' are given by 
Er=,jLy^ (z) =y^' V (3/2) = Y^ ijy) ^^ , 
(4-7r) 



a\r)=p.,{z)-A,{^)=l- 



(4.50 



4.2. Estimation 

The Rayleigh distribution is completely specified 
if the parameter y is known. If Zi, z-i, . . ., ^at is 
a sample of A^ independent observations from this 
distribution, the likelihood function (the joint prob- 
ability density of ^i, . . ., z^ consichMcHl as a 
function of 7) is given by 



where 



7 7 



1 ^ 



(4.6) 



(4.7) 



This shows that c is suflicient for 7. Since Ec^y, 
c is unbiased. We know that if an unbiased sufficient 
estunator exists, it is automatically the most eflH.cient 
estimate, also the jnaxhnum likelihood estimate. 
The variance of c is, of course, 



var c= 



11 

N' 



(4.8) 



If /. is any othei- unbiased estimate of 7, then var 
c<var t. 

Since 2Zi/y is a x^ variate with 2 degrees of freedom, 
and Zi, . . ., z^ are independent 

2Nc/y=2 S ^z-/7, 

i=l 

is a x^ variate with 2N degrees of fiecMlom. Thus 
the distribution of c is 

dP{c)= AT-p/xTx g ^ c^-Wc, if c>0; =0, otherwise. 



Vr(A^) 



(4.9) 



To set up confidence lunits on 7 with a confidence 
coefficient I — a, we determine from x^ tables two 
numbers xi and xh corresponding to 2A^ degrees 
of freedom, such that 



Then 



Pr(x'<xf)=cx/2,Pr(x'<X^ = l-a/2, 



xi<2Nc/y<x:i, 

has probability I — a and 



2Nc 



<7< 



2Nc 



(4.10) 



Xi — xt 

will be 100(1 — 0;) percent confidence Ihnits for 7. 

Example 1. Let iV=3, 1 — a = 0.98. From x^ 
tables we find, for 6 degrees of freedom, Pr 
(x'<0.872)=0.01, Pr(x'<16.812)=0.99. Hence 



T^c<y< 



6 



or 



16.812 -'- 0.872 
0.357c<7< 6.88c 
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are 98 percent confidence limits for 7. Similarly 

0.476c<T<3.67c 

are 90 percent confidence limits. 

If iV>15 so that the nmnber of degrees of freedom 

/:=2iV>30, 



2X 



-V27=T-7lf- 



'^Jm-l 



is approximately a standard normal variate [Cramer, 
1951, p. 251]. If X is a number such that 



1 r -'- 

\ Ztt J -X 



then 



4Nc 



<y<- 



4Nc 



(^l4N-l+xy- (-yj'AN-l-x) 



(4.11) 



(4.12) 



are 100 (1 — a) percent confidence limits for 7. When 
A^is large, say A7'>100, (4.12) is approximated by 



Nc 



'^JN+: 



<y< 



Nc 



■yjN—X 



(4.13) 



Example 2. Let A^=100, l-a:=0.95. We find 
0^=1.96. From (4.12) 

0.831c<7< 1.23c. 

If, on the other hand, we use (4.13), we get 

0.836c<7< 1.24c. 

Since c is the sufficient statistic for the Rayleigh 
distribution, if it is known that the sample has come 
from a Rayleigh distribution, we simply evaluate 
c][and estimate the entire distribution b}^ 



P(0) = 1 — g ~c^ 0<^<oo. 



(4.14) 



The usual practice of estimating the population 
percentiles by the sample percentiles is utterly un- 
justified. The distribution percentile, f^, should be 
estimated by 



S(?=— cln (1 — g), 



(4.15) 



which is an unbiased and minunum variance esti- 
mator of ^q. Thus, for example, the median, fg, 
will be estimated by 

s.5=-c hi 2=0.69315c, 

which has the variance 

0.4807^ 



var 0.5= (hi 2)- var c= 



N 



On the other hand, the asymptotic variance of the 
sample median, ^"5, is given by 



var 0.5 



1 



Thus 



4Ar/(f.5) N 
var 3.5 



var s. 5 



^0.480, 



i.e., the sample median is only 48 percent as efficient 
as 3.5. 

To take another example: the first decile, f.i, is 
most efficiently estimated by 

^i=--cln (.90)=0.10536c, 

which is unbiased and has variance 

var 3. 1=0.0111 — 

The sample 1st decile, z[i is biased and has 
asymptotic variance 



, 0.1X0.9 72 

^ar^.i=^^y-27 



so that 



var 3.1 



var 



<^.i 



Nf{^.,) m' 



^9X0.0111^0.1, 



i.e., the efficiency of z[i is only 10 percent. 



4.2.1. Estimation cf the Amplitude Characteristics 

We know that if sufficient statistics are available 
for a distribution, its characteristics are most effi- 
ciently estimated in terms of these statistics. Thus 
the characteristics of the amplitude distribution are 
estimated in terms of c. 

For example to estimate the mean amplitude, 
(x7)^^^/2, we need only to find a function of c which 
is an unbiased estimator of this quantity. 

Now 



Ec^=E{y+c-^y^ 

c— 7 1 {c — ^Y 



=te\i+\ 



2 7 



•] 



'[-^} 



Thus, to the order N-\ [l + l/(87V)]c^/2ig ^^ unbiased 
estimate of 7^^^ The mean of r is then estmiated by 



m = 






mP^'- 



(4.16) 
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Since 



var c^=£'c— (£'c^)2 
=7— ?( 



1- 



8A7 ^4iV' 



var 7n^- 



7r7 



On the other hand, let rt = z\, aud 



(4.17) 





- 1 #v 


We have 






- var r,- 7(4— tt) 


Thus 






varm_ x 



var r 4(4— tt)" 

i.e., the sample nu^an, T, is oiilx' 91 percent as efficient 
as the unbiased estiinalcir based on the sufficient 
statistic c. 

Similar considerations apply to estimating other 
characteristics of the amplitude distribution, and 
hence there is no need for us to go into derails. 

4.2.2. Estimation When the Data Are Truncated 

Sometimes, due to the limitations of the recording 
instrument, observations above a certain level x are 
not recorded. However, the prior information is 
there that the unrecorded observations will be greater 
than X, and the recorded ones will be less than x. 
This prior information together with the actual values 
of the recorded observations should be utilized to 
obtain the most effi,cient estim.ator of 7. 

In a sam.ple oi N independent observations let N—k 
values exceed x, and let the remaining k observa- 
tions be ^1, 2-2, . . ., 2]c. Obviously, k itself is a 
random variable with possible values 0, 1, 2, . . ., 
N, Also the probability density function of Zi is now 
conditioned upon the knowledge that 0<Zi<x, 

The conditional probability density function is 

1 



j)(z\{)<z<x)^-~e- 



-""',^<z<x, 
The probability distribution of k is 

^ )v\i-vY-\k=QA, 

Hence the likelihood function is 



wherep = l— 6 "> 

PriK=k) = 



,,N. 



L{y)- 



-Nx 



~,iik=0, 



0<z„, . .,z,<x,k=l,2,. . .,N. 



(4.18) 



Obviously, when k = no point esthnate of 7 is 
possible, since the only information available is that 
all observations are greater than x. In this case, 

which has probability 6 y , we can oidy infer that 
x<iy<^ °° . When ky^O, the maximum likelihood esti- 
mate of 7 is 

{N-k)x+'^Zi 
ci= ^ '^' > (4.19) 

The properties of Ci are not as easy to establish as 
those of c in section 4.2. Using the well known 
properties of the maximum likelihood estimators, 
however, we can say that for large A^, Ci will be 
normally distributed with mean Eci, and variance, 
var Ci, where 

^ , x(l—p) 



var Ci^ 



27^ 



Np' 
{4p- 



-3)r 



Np Np^ 



,p=l-e y. (4.20) 



4.3. Is the Variate Rayleigh Distributed? 

Let Zi, Z2, . . ., 2v be independent observations on 
a non-negative variate z. We wish to test whether 
z is Rayleigh distributed. Proceed in the following 
manner. 

(Choose an integer m, not too small (at least 5), 
such that N/'tn is greater than 5, preferably greater 



N 

than 10. Calculate c=^Zi/N. 

i=\ 



Deterj 



nnu^ luini- 



bers Xi, 



Xm-i, such that 



P(xd = l- e'=~,i=l,2, 



m- 



i.e. 



Xi^c In ( : \ i=l, 2, . . ., m— 1. 

\m—ij ' * 

The range (0,oo) of z is then divided into m nonover- 
lapping intervals /i = (0,Xi), l2={xi,Xy), . . ., /^= 
(x^_i, 00)^ such that the expected number of obser- 
vations, under the hypothesis F{z), in each interval 
is NIm. 

Let the number of actual observations falhiig in 
these intervals be /i,/2, . . .,fni- Then 



xHS 



N/m 



7Y) ^^* 



■N (4.21) 



is approximately a x^ variate with m— 2 degrees of 
freedom [Cramer, 1951, pp. 424-426]. Pre-assign 
some critical probability level a:(=0.05, or 0.01, say) 
and let xi b^"^ the number such that 

Prix'>xi)=a, 

If the observed value Xo<Xaj accept the hypothe- 
sis that 2 is a Rayleigh variate, and estimate its dis- 
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x§>x^ 



reject the 



ti'ibution function by P{z). If 
hypothesis of Rayleigh distribution. 

If the sample is truncated, i.e., if the maxin).um 
value which could be recorded is x, we estimate y by 



es are 



c^=\ {N—k)x+^Zi\k, where A^— ^ valu 

greater than x and unrecorded. We proceed as 
above, changing the last interval (or intervals), if 
necessary, to read (:c,oo) and determining the 
expected number of observations in this interval by 






fo=oh- 



Ne-'^l^i] and calculating xl=^ 

served, /e= expected frequency. 

Example 3. The following is a systematic sample 
of 80 observations (read at 5 sec intervals) of re- 
ceived field intensity in (microvolts) ^ Tlie follow- 
ing are the transmission parameters: 

Frequency: 400 Mc/s 

Transmitter power: 2 kw 

Angle : zero 

Path length : SOOmilliradians (165 miles) 

Antennas: 60 ft diam parabolic re- 

flectors having a gain of 
38 db over an isotropic 
radiator 

The measured hourly me- 
dian basic transmission 
loss ^209 db 
0425 a.m., July 27, 1963 



Loss: 



Time: 



Observed values of received power in (mv)^ (read 
left to right) : 

0.05, 0.76, 0.32, 0.96, 0.63, 0.09, 

0.26, 0.10, 0.95, 0.01, 0.50, 1.26, 

0.01, 0.16, 0.56, 3.16, 1.27, 2.24, 

0.28, 0.21, 0.35, 0.20, 0.39, 0.89, 

1.01, 0.49, 0.90, 1.90, 1,42, 1.56, 

2.40, 2.24, 0.80, 0.56, 1.45, 0.18, 

0.18, 1.31, 0.64, 1.95, 0.48, 0.55, 

0.71, 0.48, 0.40, 0.06, 0.79, 1.01, 

0.16, 0.01, 0.06, 0.03, 0.01. 



0.20, 
0.18, 
1.99, 
1.00, 
1.24, 
1.32, 
0.02, 
0.44, 



0.71, 
0.25, 
0.32, 
0.81, 
0.08, 
1.20, 
0.28, 
0.28, 
0.51, 0.70, 0.14, 

We have N= 



0.06, 
0.45, 
0.51, 
1.29, 
0.98, 
1.59, 
0.81, 
0.07, 



80, c = 0.71. Take m = 8, a=0.05. 



For 6 degrees 
is 12.592. 



of freedom the critical value of x^ 



Calculate :ri = --0.71 In 



^-l, 2, 



7. 



The resulting class intervals and the observed 
frequencies are as follows. 



Class intervals 


fo 


-0. 095 


13 


0. 095- . 204 


9 


. 204- . 334 


8 


. 334- . 492 


8 


. 492- . 696 


8 


. 696- . 984 


13 


. 984-1. 476 


12 


1. 476- oo 


9 



The expected frequency in each interval is 10. We 
have 

xH^S/'-80 = 3.G<12.592. 

The data are consistent witli the hypothesis that 
z is Rayleigh distributed. In fact, the probability 
that such a sample or worse, as measured by x^, 
comes from, a Rayleigh distribution is inore tlian 
70 percent. The entire (hstrihution is estiirated by 



F{z) = l-e 



0<2<oo. 



4.4. Two Independent Samples 

Sometimes we wish to test the hypothesis that 
two samples, which are known to be ivoin Rayleigh 
distributions, are from the same distribution, i.e., 
7i=72. For example, to test the assumption of 
stationarity of a Rayleigh process, we may take 
two samples from the record (sufficiently far apart 
in time to insure independence), and test the hypoth- 
esis that 7i=72. Or, we may be recording a 
signal, not continuously, but after breaks in be- 
tween recordings, and may wish to test whether 
the characteristics of the received signal have 
changed from record to record. 

Let it be known that two independent samples, 
^1, 02, • • • , Zn^ and 2i, 09, . . . , z\f, are from 
Rayleigh distributions. To test tlie hA'pothesis 
that 7i=72, it is sufficient to note tliat 



F{2N, 2M) . 



N 

~~ 2N 



2M_ 

M 

1=1 



'C2 



is a Fisher-Snedecor F variate with the indicated 
degrees of freedom. Here, Ci and C2 denote the 
sample means of the first and second sample re- 
spectively. We assume that Ci > C2] if not we 
simply invert the ratio and interchange the degrees 
of freedom. Pre-assign a significance level a, and 
test for the significance of the calculated F. 

Example 4. Consider the data in example 3. Let 
us test the hypothesis that the first 25 observations 
are from the same Rayleigh distribution as the last 
25 observations. Let the significance level be 
taken as 0.05. From F tables the 5 percent significant 
value, for 50 and 50 degrees of freedom, is 1.60. 
Denoting by Ci and C2 the mean values of the first 
25 and the last 25 observations, respectively, we have 

Ci=0.578, C2=0.482. 
Hence, the calculated value of Fis 

F ^^^ 1 90 
^^=482^^-^^^ 

which is less than 1.60. Hence we accept the 
hypothesis that 71=72. 
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5. Distribution of the Ratio of Two Rayleigh 
Variates 

Manx" pi-obleins in radio engineering require the 
knowledge of the distribution of the ratio of two 
independent Rayleigh variates, such as the signal 
power and the noise power, or the powers of two 
interfering signals. 

Let Zi and Z2 be two independent Rayleigh variates 
with parameters 71 and 72. Let u=Zi/z2, and 
€=71/72. Since 2^1/71 and 2^2/72 are independent 
X^ variates each with 2 degrees of freedom, 



^i/Ti 






^2/72 



€32 



is a Fislier-Snedecor /" variate with 2 and 2 degrees 
of freedom [Cramer, 1951, pp. 241-242]. Hence the 
distribution of V;=eUi is given by 



dP(u)-- 



edu 



{e+uy 

l-P{u)=Pr{zJz2>u) 



0<u< 00; 



e+u 



(5.1) 



(5.2) 



Example 5. Let the average signal ])()wer be 20 db 
above the average noise power so that e = 100. Let 
satisfactory service require signal powder 10 db above 
the noise power. To calculate the probability of 
satisfactory service we set € = 100, 1^ = 10 in (5.2) and 
obtain Pr(si/s2>10)=0.91. 

The distribution of u involves only one parameter 
€=71/72, which is the median of the distribution. 
We note that the mean, variance, and higher 
moments of u do not exist. 

Let Ui, . . .^ Un be independent observations on 
u. The logarithm of tlie likelihood function is 



lnL(0=iVln€-2X;in {e+u^, 



(5.3) 



which sliows that no sufficient estimator exists for e. 
The maxhnum likelihood esthnator, e, of € is the 
solution of the equation 



blni 



=0, 



i.e., of 



iV_^_2__ 
e i=ie+Ui 



(5.4) 



This equation is quite difficult to solve if A^>3. 
We, therefore, find some less efficient but easily 
available esthuator. The Cramer-Rao greatest 
lower bound for the variance of any unbiased 
esthnator is found to be 



-H"^)']-- 



N' 



(5.5) 



which is also the asymptotic variance of e. The 
asymptotic variance of the sample inedian, u^^, 
is given by 



var u^5^ 



4Np'(e)~N' 



(5.6) 



so that the efficiency of the sample median is 75 
percent. Thus, for a quick esthnate of e, the median 
may be used. In any case the sample mean should 
not he taken as an estimate of e. 



6. A Rayleigh Process Derived From a 
Normal Process 

Referring back to the tli(M)retical framework of 

section 3, w^e assume now tliat the resultant vector r 

is a function of tune, i.e., r{t) is a random vector 
process. The components (x, y) will now be written 
{x{t), y(t)) and it will be assumed that x(t) and y{t) 
are independent stationary normal processes with 
the same variance a^ and the same autocorrelation 
function p(s). That is to say, for all t and s 

ExHt) = Ef(t)=a', 

Ex(t)x {t+s)=Ey {t)y {t+s) = a'p(s) . (6.1) 

z(t)=x'(t)+f(t), y=2<y\ (6.2) 

luMi, for all / and s^ 

Ez(t)=y, 

Ez{t)z{t+s)=y'{l+p'(s)). (6.3) 



Set 



Wj'iting a(s) for the autocorrelation of tlie z{t) 
process, we have 



,, cov(zit),z(t+s)) , 

a{s)= ^ ^p {s). 

7 



(6.4) 



The process is completely determnuMl if 7 and 
a(s) = p^(s) are known. 

As before we use the sample mean now giviMi by 

c=^j%{t)dt, (6,5) 

to estimate 7. c is unbiased and has variance 

varc=^J (l--^a{s)ds. (6.6) 

We note that in the case of a discrete sample 

z(l), z(2), . . ., ziN), (6.6) will read 



y2 272A^-1 



iO-^)« 



(S). 



(6.6') 
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The exact distribution of c is unknown. However, 
let A^' be defined by the equation 

var c=77A^', i.e., A^'=T7var c, (6.7) 

where var c is given by (6.6) or (6.60- N^ will be 
called equivalent random sample size for estimating y. 
The results of sections 4.2-4.4, then, will approxi- 
mately hold with A^ and M replaced by equivalent 
random sample sizes A^' and M\ N\ in general, 
will not be an integer, but this presents no difficulty 
either in theory or in calculations. If [2A^'] is the 
largest integer in 2N\ we will sinaply interpolate 
between the percentiles of x^ for [2N'] and [2iV'] + l 
degrees of freedom, to obtain the percentile points 
of 2N'c/y. Similarly, for F(2N\2M') of section 4.4. 



Example 6. Let a{s)^a{ — s)-- 
Also, let p=e~^, so that a(s)=p\ 



=:6-^^ 8>0, iu>0. 
From (6.6) 



var c 



=r;[-^> 



V 
mT 



if /x Tis large. Hence, the equivalent random sample 
sizeN'^fjLT/2. 

In case of a discrete sample, we obtain from (6.6^ 



var 



'4m> 



so that N'^N{l-p)/(l+p). 

For example, if e7'^=:p = 0.5, m = 0.69315, 2N'-- 
0.6931 5 T in the continuous sample case, and 2A^' = 
2A^/3 in the discrete sample case. 



Let 



Then 



C(s)^^J^^~' z(t)z(t+s)dt. 



a(s)=- 



(6.8) 



(6.9) 



is a consistent estimate of a(.s'). 

From the joint distribution of z(t) and ^(^+6), 
obtained from Siddiqui's [1961] equations (5.8) and 
(5.9) by setting n=2, we also find the following: 

E(z(t+s)\z{t))=a{s)z{t)+y(l-a(s)), (6.10) 

so that the regression of z(t-]-s) on z(t) is linear. 
However, 

V^v(z{t+s)\zit))=y\l-ais)y+2ya(s){l-ais))z(t), 

(6.11) 
is not independent of z{t). 
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